Search for: All records

Abstract The symmetry breaking in octahedral silsesquioxane and its germanium analogues (Si₈O₁₂H₈and Ge₈O₁₂H₈) has been investigated using the M06-2X/6-31++G(3df, 3pd) method and group theory. Both structures undergo$${O}_{h}\downarrow {T}_{h}$$symmetry breaking, characterized by pseudo-Jahn−Teller stabilization energies of 0.22 kcal/mol for Si-POSS and 9.82 kcal/mol for Ge-POSS. Under the influence of the pseudo-Jahn–Teller effect, the distortion vector involves the vibrational a_2gmode with imaginary frequency. The distortion forces in O_h-POSS are predominantly localized on the oxygen atoms and driven by the coupling between the lowest unoccupied molecular orbital (a_1g) and the highest occupied molecular orbital (a_2g). The symmetry breaking is attributed to a pseudo-Jahn–Teller mechanism of type (a_2gx a_1g) = a_2g. The symmetrical substitution of oxygen atoms by X (where X = C, N, P) results in viable T_h-Si₈X₁₂H₈and T_h-Ge₈X₁₂H₈structures. The observed pseudo-Jahn–Teller distortion and substitutional symmetry breaking caused by X indicates a consistent electronic relaxation mechanism, characterized by the formation of C=C, N=N and P=P bonds on the POSS cubic faces, which serves as hallmarks of stability. Additionally, we find that the volume of substituted T_h-symmetrical POSS is sufficiently large to accommodate small ions. 

ABSTRACT The energies and geometries of the lowest lying singlet and triplet states of the four diradicals formed by removing two H atoms from thiophene have been characterized. We utilized the highly correlated, multireference methods configuration interaction with single and double excitations with and without the Pople correction for size‐extensivity (MR‐CISD+Q and MR‐CISD) and averaged quadratic coupled cluster theory (MR‐AQCC). CAS (8,7) and CAS (10,8) active spaces involving σ, σ*, π, and π* orbitals were employed along with the cc‐pVDZ and cc‐pVTZ basis sets. The larger active space included the two electrons in the nonbonding sp²hybrid orbital on sulfur. We find that all didehydro isomers exist as planar, stable ground state singlets. The singlet‐triplet (S‐T) adiabatic gaps range from 15 to 25 kcal/mol while the vertical splittings are 21–35 kcal/mol. The 2,3 isomer has the lowest absolute ground state singlet energy and the largest adiabatic and vertical S‐T splitting. The ground states of the 2,3‐, and 2,5‐didehydrothiophene isomers are predicted to exhibit the smallest and largest diradical character, respectively, based on their electronic structures, spin densities and bonding analysis. To our knowledge, no experimental excitation energies of any of the didehydrothiophene isomers are available, and our computed MR‐AQCC/cc‐pVTZ data are believed to be among the most accurate computed results. This extensive study shows a competitive performance between MR‐AQCC and MR‐CISD+Q. 

Heteroaromatic species are commonly found in complex gaseous mixtures, from tobacco smoke to petroleum and asphaltene combustion products. At high temperatures, C–H bond rupture produces various dehydro radical isomers. We have used the spin–flip formulation of equation-of-motion coupled cluster theory with single and double substitutions (EOM-SF-CCSD) to characterize the energies and wave functions of the lowest lying singlet and triplet states of the diradical (2,3), (2,4), (2,5), and (3,4) di-dehydro isomers of pyrrole, furan, and thiophene. In all cases, these diradicals are minima on the broken-symmetry ωB97X-D/cc-pVDZ potential energy surface. In most cases, the diradical geometries distort to enhance through-space or through-bond coupling in the singlet states and to avoid Coulombic or exchange repulsion in the triplet states. EOM-SF-CCSD results indicate that all diradical isomers are two-configurational, closed shell singlet states. The only exceptions to this are for (2,3) and (2,4) thiophene and (2,3) pyrrole, which each contain more than two configurations. In all cases, the leading term in the multiconfigurational diradical wave function doubly occupies the symmetric radical σ orbital, indicative of either through-space or 1,3 through-bond coupling. We utilized the nucleus-independent chemical shift (NICS) approach to qualitatively assess aromaticity and find that this property varies and may be related to the energetic splittings in these diradical isomers. 

Abstract Oligomeric models of linear ladder silanes, siloxanes and siloxazanes with seven repeat units consisting of four-, six-, or eight-membered rings were designed and their conformations in chloroform were explored. The Low Mode–Monte Carlo conformational method was used to explore oligomeric flexibility on the OPLS-2005/GBSA(CHCl₃) potential energy surface to obtain a set of low energy structures for each oligomer. These structures were then optimized using B3LYP/6-31G*/SCRF-PBF(CHCl₃) calculations. The results indicate complex conformational dynamics with mostly non-planar, curved structures. Electron delocalization from the lone pair of electrons on N or O into empty 3d orbitals on Si was not observed. 

The COVID-19 pandemic was declared due to the spread of the novel coronavirus, SARS-CoV-2. Viral infection is caused by the interaction between the SARS-CoV-2 receptor binding domain (RBD) and the human ACE2 receptor (hACE2). Previous computational studies have identified repurposed small molecules that target the RBD, but very few have screened drugs in the RBD–hACE2 interface. When studies focus solely on the binding affinity between the drug and the RBD, they ignore the effect of hACE2, resulting in an incomplete analysis. We screened ACE inhibitors and previously identified SARS-CoV-2 inhibitors for binding to the RBD—hACE2 interface, and then conducted 500 ns of unrestrained molecular dynamics (MD) simulations of fosinopril, fosinoprilat, lisinopril, emodin, diquafosol, and physcion bound to the interface to assess the binding characteristics of these ligands. Based on MM-GBSA analysis, all six ligands bind favorably in the interface and inhibit the RBD–hACE2 interaction. However, when we repeat our simulation by first binding the drug to the RBD before interacting with hACE2, we find that fosinopril, fosinoprilat, and lisinopril result in a strongly interacting trimeric complex (RBD-drug-hACE2). Hydrogen bonding and pairwise decomposition analyses further suggest that fosinopril is the best RBD inhibitor. However, when lisinopril is bound, it stabilizes the trimeric complex and, therefore, is not an ideal potential drug candidate. Overall, these results reveal important atomistic interactions critical to the binding of the RBD to hACE2 and highlight the significance of including all protein partners in the evaluation of a potential drug candidate. 

Undergraduate research and education are critical parts of the scientific ecosystem. Professors at primarily undergraduate institutions ({PUIs}) bear responsibility for molding young scientists while carrying out high quality research. In this Voices piece, researchers offer inspiring replies to the following prompt: in your role as a professor, what are the challenges and opportunities when conducting research at a {PUI}?},